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ABSTRACT 

Using the semi-analytic method proposed by Esmailzadeh et al. (1991) we 
calculate the abundances of the light elements produced during primordial nu- 



Ph' cleosynthesis assuming that the gauge coupling constants of the fundamental in- 

teractions may vary. We analyze the dependence of the nucleon masses, nuclear 
binding energies and cross sections involved in the calculation of the abundances 
with the fundamental constants assuming the chiral limit of QCD. The abun- 
dances of light elements as a function of the fundamental constants are obtained. 
Finally, using the observational data of D, 3 He, 4 He and 7 Li we estimate con- 
straints on the variation of the fundamental constants between the primordial 
nucleosynthesis and the present. All observational abundances and the WMAP 
estimate of the baryon density, can be fitted to the theoretical predictions with 
varying coupling constants. The possible systematic errors in the observational 
data, precludes from stronger conclusions. 

Subject headings: primordial nucleosynthesis, varying fundamental constants, 
cosmology 



1 Facultad de Ciencias Astronomicas y Geofisicas. Universidad Nacional de La Plata. Paseo del Bosque 
S/N 1900 La Plata, Argentina 

2 Departamento de Fi'sica, FCEyN, Universidad de Buenos Aires, Ciudad Universitaria - Pab. 1, 1428 
Buenos Aires, Argentina 



- 2 - 



1. Introduction 

Big Bang Nucleosynthesis (BBN) is one of the most important tools to study the early 
universe. The model is simple and has only one free parameter, the density of baryonic 
matter, which can be determined by comparison between theoretical calculations and ob- 
servations of the abundances of the light elements. On the other hand, data on cosmic 
microwave background (CMB) provide an alternative, independent method for determining 
Vt B h 2 (Spergel et al. 2003). Recently, the concordance between both methods has been inves- 
tigated by many authors (Cyburt et al. 2003; Romano et al. 2003; Cuoco et al. 2004; Cyburt 
2004; Coc et al. 2004a,b). From the WMAP baryon density , the predicted abundances are 
highly consistent with the observed D but not with 4 He and 7 Li. They are produced more 
than observed. Such discrepancy is usually ascribed to non reported systematic errors in the 
observations of 4 He and 7 Li. Indeed, more realistic determinations of the 4 He uncertainty 
implies a baryon density in line with the WMAP estimate (Cyburt 2004; Olive and Skillman 
2004). On the other hand, Richard et al. (2005) have pointed out that a better understanding 
of turbulent transport in the radiative zones of the stars is needed for a better determination 
of the 7 Li abundance. However, if the systematic errors of 4 He and 7 Li are correctly esti- 
mated, we may have insight into new physics beyond the minimal BBN model, for example: 
new neutron lifetime (Mathews et al. 2005), super WIMP scenario (Feng et al. 2003), lepton 
asymmetry (Ichikawa et al. 2004) and varying constants (Bergstrom et al. 1999; Nollet and 
Lopez 2002; Ichikawa and Kawasaki 2002, 2004). Therefore, BBN is not only one of the most 
important tests of the Big Bang theory, but it is also useful to obtain stringent constraints 
on deviations from standard cosmology and on alternative theories to the Standard Model 
of fundamental interactions (SM). 

Among these theories, there are some in which the gauge coupling constants may vary 
over cosmological time scales like string derived field theories (Wu and Wang 1986; Maeda 
1988; Barr and Mohapatra 1988; Damour and Polyakov 1994; Damour et al. 2002a,b), related 
brane- world theories (Youm 2001a,b; Palma et al. 2003; Brax et al. 2003), and (related or 
not) Kaluza-Klein theories (Kaluza 1921; Klein 1926; Weinberg 1983; Gleiser and Taylor 
1985; Overduin and Wesson 1997). On the other hand, a theoretical framework in which 
only the fine structure constant varies was developed by Bekenstein (1982) and improved 
by Barrow et al. (2002). This model was generalized in order to study the time variation 
of the strong coupling constant (Chamoun et al. 2001). Different versions of the theories 
mentioned above predict different time behaviors of the gauge coupling constants. Thus, 
bounds obtained from astronomical and geophysical data are an important tool to test the 
validity of these theories. 

The experimental research can be grouped into astronomical and local methods. The 
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latter ones include geophysical methods such as the natural nuclear reactor that operated 
about 1.8 10 9 years ago in Oklo, Gabon (Damour and Dyson 1996; Fujii et al. 2000, 2002), 
the analysis of natural long-lived f3 decayers in geological minerals and meteorites (Dyson 
1966; Sisterna and Vucetich 1990; Smolliar 1996) and laboratory measurements such as 
comparisons of rates between clocks with different atomic number (Prestage et al. 1995; 
Sortais et al. 2000; Marion et al. 2003). The astronomical methods are based mainly in the 
analysis of spectra form high-redshift quasar absorption systems (Cowie and Songaila 1995; 
Varshalovich et al. 1996; Webb et al. 1999, 2001; Murphy et al. 2001a,b; Levshakov et al. 
2002; Ivanchik et al. 2002; Murphy et al. 2003; Ivanchik et al. 2003; Bahcall et al. 2004). 
Although, most of the previous mentioned experimental data gave null results, evidence of 
time variation of the fine structure constant was reported recently from high-redshift quasar 
absorption systems (Webb et al. 1999, 2001; Murphy et al. 2001a,b, 2003; Ivanchik et al. 

2003) . However, other recent independent analysis of similar data (Martinez Fiorenzano 
et al. 2003; Quast et al. 2004; Bahcall et al. 2004; Srianand et al. 2004) found no variation. 
On the other hand, measurements of molecular hydrogen (Ivanchik et al. 2002, 2003) reported 
a variation of the proton to electron mass /i = ^ 

The time variation of the gauge coupling constants in the early universe can be con- 
strained using data from the Cosmic Microwave Background (CMB) (Battye et al. 2001; 
Avelino et al. 2000; Martins et al. 2002; Rocha et al. 2003) and the primordial abundances of 
light elements (Bergstrom et al. 1999; Nollet and Lopez 2002; Ichikawa and Kawasaki 2002, 

2004) . 

The prediction of the light elements abundances ( 4 He, D, 7 Li) produced during the 
first minutes of the universe can be calculated using numerical (Wagoner 1973; Kawano) 
and analytical (Esmailzadeh et al. 1991; Mukhanov 2003) methods. Ichikawa and Kawasaki 
(2002) modified the public code in order to analyze the BBN scenario with varying gauge 
coupling constants. They considered a theoretical model taken from string theory where the 
variation of the coupling constant is related to the expectation values of the dilaton field and 
compared with observational data. In consequence, the results they obtained are restricted to 
the validity of this model. Furthermore, numerical calculations of the theoretical abundances 
of the light elements allowing only a variation of the fine structure constant were performed 
by different authors (Bergstrom et al. 1999; Nollet and Lopez 2002; Ichikawa and Kawasaki 
2004). On the other hand, an analytical study of 4 He abundance including variation of 
the gauge coupling constants was performed by Miiller et al. (2004). Moreover, the change 
in the abundance of 4 He due to variable mass in 5 dimensional theories was analyzed by 
Anchordoqui et al. (1996). Finally, the effect of considering non extensive thermostatistics 
has been analyzed by various authors (Torres et al. 1997; Pessah et al. 2001; Pessah and 
Torres 2001). 
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In this work, we follow the semi-analytical method proposed by Esmailzadeh et al. (1991) 
to study the effect of a possible variation of the values of the three gauge coupling constants 
of the Standard Model of Particles Interactions (SM) between primordial nucleosynthesis and 
the present. Even though, the semi-analytical method gives results one order of magnitude 
less accurate that the calculations performed with the numerical code, it is very useful to find 
out the dependence of the abundances and temperatures with the fundamental constants, 
which is one of the principal aims of this work. 

We will not assume any of the theoretical models for varying constants mentioned above. 
Motivated by theoretical predictions and observational data, we will study the formation of 
the light elements in the early universe assuming that the values of the gauge coupling con- 
stants of the fundamental interactions (electromagnetic, strong and weak) may be different 
from their actual value. Thus, our approach is a phenomenological one and our results will 
be model independent. Furthermore, we assume the chiral limit of QCD to analyze the de- 
pendence of nucleon masses, binding energies and cross sections with the strong interaction 
coupling constants. The gauge coupling constants of U(l), SU(2) and SU(3), namely, a±, 
a 2 and a 3 are related with the fine structure constant a, the QCD energy scale Aqcd and 
the Fermi coupling constant Gp through the following equations: 

a' 1 (E) = ^r 1 (E) + a^ 1 (E) (1) 



Aqcd = E exp 



2 

2tt 



y % 1 (E) 



(2) 



71 a 2 (Mz) ( „, 

Gf = -^uf (3) 

where E refers to the energy scale and M z refers to the boson Z mass. Actually, we will 
study the dependence of the different physical quantities involved in the calculation of the 
primordial abundances with a, Aqcd and Gf- 

Almost all of the observational and experimental data are consistent with no variation 
of the constants (Landau and Vucetich 2002). Moreover, the reported variations (Murphy 
et al. 2003; Ivanchik et al. 2003) are very small (^ i ~ 10~ 5 ). Therefore, in order to find 
out the dependences of relevant physical quantities with a, Aqcd and Gf we will perform a 
Taylor expansion to first order in each case as follows: 

QQ 

AO = — /„ioda« /.today ^today. Aft 
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where Q refers to the physical quantities involved in the nucleosynthesis calculation such as 
nucleon and nucleus masses, nuclear binding energies, cross sections and abundances of the 
elements. 

In the standard picture, the only free parameter of the nucleosynthesis calculation is 
the density of baryonic matter Qsh 2 . This quantity has been determined with a great 
accuracy with data from the CMB provided by WMAP (Spergel et al. 2003). On the other 
hand, the baryon density can also be estimated using data provided by galaxy surveys 
(SDSS, 2dF) and x-ray satellites (Chandra, XMM- Newton, ROSAT, ASCA). In appendix 
A we combine different data to obtain an estimation of fish 2 independent of the WMAP 
estimate. Therefore, we shall approach to the problem studying the dependences of all 
physical quantities and abundances with both the fundamental constants and the baryon 
density. Thus, we will obtain the uncertainties of the abundances of the light elements as 
function of the variations of the fundamental constants with respect to their actual value 
and as function of the variation of Qsh 2 with respect to the WMAP estimate (Spergel et al. 
2003). On the other hand, we will also compare the predicted theoretical expressions for 
the abundances with observational data and include independent estimates of the baryon 
density in the analysis (see section 5). 

Furthermore, in section 2, we shall calculate the dependence of the nucleon masses 
and binding energies with the fundamental constants, and in section 3, the corresponding 
dependence of the relevant scattering cross sections. We have carried this calculations in 
some detail, since there are several subtle points in these dependences that will be clearly 
exhibited in the final results. In section 4 we apply the semi-analytical method proposed by 
Esmailzadeh et al. (1991) to calculate the abundances of the light elements and their de- 
pendence with the fundamental constants. In section 5 we briefly describe the observational 
data and the results of comparing them with the theoretical predictions calculated in this 
work. We also discuss our conclusions. 



2. Masses and binding energies of light elements 

In this section we analyze the dependence of the nucleon masses, nuclear binding energies 
and nuclei masses of the light elements with the fundamental constants a and Aq C d- The 
weak interaction contribution is too small to produce any observable consequences (Haugan 
and Will 1976; Chamoun and Vucetich 2002). 

The dependence of the hadronic masses and nuclear binding energies with the QCD 
coupling constant a 3 or the QCD scale parameter Aq C d depends on the model of hadronic 
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interactions considered. However, if we assume that the quark masses are null, an assumption 
which is called in the literature as chiral limit, there is a only a single parameter in the 
theory, namely the QCD scale parameter Aqcd- Even though great efforts (Beane and 
Savage 2003b, a; Epelbaum et al. 2003; Flambaum and Shuryak 2002, 2003; Dmitriev and 
Flambaum 2003; Dmitriev et al. 2004; Olive et al. 2002) have been done in order to analyze 
the dependence of nucleon masses and binding energies with Aqcd beyond the chiral limit, 
this task is not trivial and highly model dependent. 

On the other hand, from simple dimensional analysis (Stevenson 1981), it follows that 
in a theory with only one relevant parameter all static observables with dimension of mass 
must be proportional to this parameter, which in our case is Aqcd- More precisely, any 
quantity a with units of E D (where E means energy) must satisfy an equation of the form: 



° = A QCDf 



Q 



A 



QCD 



(5) 



where Q is a quantity specifying the energy scale. Furthermore, for static quantities such as 
nucleon masses the previous equation takes the form: 



a qcdI 



Aqcd 



(6) 



since the only scale parameter is a itself. The solution of equation 6 reads: 

° = KcdX (7) 

where X is a dimensionless numerical constant. In such way, all low-energy static quantities 
will satisfy an equation of the form 7. Moreover, all nucleon masses and energies will have 
a linear dependence: 

m N ~ eB ~ Aqcd (8) 

and all nuclear radii will satisfy: 

R ~ A QCD (9) 

since we use units where H — c — 1 for this analysis. The chiral limit was previously 
considered by Sisterna and Vucetich (1990) studying time variation of fundamental constants 
in planetary phenomena. 

The mass of the nucleons can be written as a sum of two contributions: the electromag- 
netic contribution and the strong interaction contribution m^: 

m N = m% + m s N (10) 



-7- 



The electromagnetic contribution depends on the nuclear radius R as follows: 

Z e 2 

Therefore the electromagnetic contribution to the nucleon mass in the chiral limit has the 
following dependence with Aqcd- 

Aqcd (12) 



Cottingham (1963) used perturbation theory to calculate the electromagnetic self energy 
of a nucleon uin to first order in a: 

Ka (13) 



where K can be expressed as a function of Sachs form factors G^ iM , which can be calculated 
from measurements of electron-nucleon scattering. On the other hand, the strong interaction 
contribution to the mass in the chiral limit is proportional to Aqcd- Therefore, we can write: 

7Tl N - m N atoday A today l i4 J 

< = (15) 

QCD 

After performing a Taylor expansion to first order, as explained in section 1 and using 
equations 14 and 15, we obtain the dependence of the nucleon masses with the fundamental 
constants: 

5m N _m%5a | 5A QCD _ p 5a | 5A QCD ^ 
m N m N a A QCD a A QCD 
The values of P are shown in table 1. 

Next, we analyze the dependence of the nuclei masses with a and Aqcd- As we did for 
nucleons, we perform a Taylor expansion to first order to obtain for a nucleus of mass m x 
the following expression: 

— - = (A-Z) — — - + Z^- — ^ (17) 

m x m x m n m x m p m x e x 

In the more general case, the binding energy (e x ) can be written as a sum of two terms: the 
electromagnetic contribution (ec) and the strong interaction contribution (es) as follows: 
£ x = + e 5- However, in the cases of nuclei with only one proton (D and T), there is no 
electromagnetic interaction and therefore the electromagnetic contribution (ec) is null . On 
the other hand, the same arguments that were used to obtain equations 14 and 15 can be 
applied for the binding energy to obtain: 

5e x e c Sa 5A QC d /1oA 
— = h — [la) 

£x £x ol Aqcd 
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Table 1: Dependence of nucleon and nuclei masses with the fundamental constants: 

pSa _|_ <?Aqco 



Nucleon / Nucleus 


P(xlO" 4 ) 


m p 


6.71 


m n 


-1.38 


D 


2.67 


T 


1.32 


3 He 


1.05 


4 He 


0.66 


6 Li 


1.50 


7 Li 


1.14 


7 Be 


2.30 



Inserting this last expression in equation 17, we obtain the general expression for the depen- 
dence of a nucleus mass with a and Aqcd- 

Sm x = p 5a | SAqcd ^ 

m x a Aqcd 

The values P for different nuclei are shown in table 1. 



3. Thermonuclear reaction rates 

In this section we calculate the thermonuclear reaction rates as functions of fundamental 
constants. We also show the dependence of the reaction rates with the baryon density 
p B = VL B h 2 . Following Esmailzadeh et al. (1991) we can write the thermonuclear reaction 
rate as: 

[ij -> kl] = p B N A (av) = 0.93 x lO~ 3 tt B h 2 T*N A (av) — (20) 

seg 

where a is the cross section, v is the relative velocity, p B = 0.93 x 10" 3 f2 B /i 2 T 9 3 ^Jip is the 
density of baryonic matter, N A is Avogadro's number per gram, T 9 is the temperature in 
units of 10 9 K. 

Using a Maxwell-Boltzmann distribution in velocities, the Boltzmann averaged cross 
section, (av) can be expressed as follows: 

M=(^) 3/2 J 'e-&vv(E)<Pv (21) 
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We need to find out as a function of the relative variations of the fundamental 

constants, (^-, ^x§§§ > ^§f) anc ^ S q f, b? : ^ ne re l at ive variation of the value of the baryon den- 
sity with respect to WMAP estimate Vt B h 2 = 0.0224 (Spergel et al. 2003). The temperature 
does not depend on the values of the fundamental constants, but the final temperature of 
each stage does and therefore, we can write: 

m - kl] _ 5T[ 

[ij -> kl] W T{ (av) 1 } 

where Tg = f(a, Aq C d, Gf) f° r & U the reaction rates. On the other hand, depends 
on the fundamentals constants through the masses of the nucleons and light nuclei and the 
form factor of the reactions. In the general case, there are not analytic expressions for cr(E) 
derived from "first principles". We suggest several expressions that attempt to fit cr(E), 
according to the elements in the reactions . 



3.1. Cross sections for charged particles reactions 

The cross section for charged particles reactions is given by Fowler et al. (1967, 1975); 
Wagoner et al. (1967): 

„ S(E) 
E 

where Zj is the charge of the i particle, /x = is the reduced mass, E is the energy, 

S(E) is the form factor. The dependence of the cross sections for charged particle reactions 
have been analyzed previously (Bergstrom et al. 1999; Nollet and Lopez 2002). In particular, 
Nollet and Lopez (2002) improved the analysis and studied the form factor as a function of 
a. In this paper, we use the criteria established by these authors to analyze the dependence 
of the form factor with a. 

Next, we analyze the dependence of the form factor with Aqcd using dimensional ar- 
guments and the chiral limit. The units of the cross section are cm 2 and therefore it follows 
that in a theory with massless quarks a ~ ^qcd- ^ ne om y quantity of eq. 23 that has units 
is the factor and thus we obtain: 

S(E) ~ Aqcd (24) 

This is valid for all charged particle reactions. The exact dependence of the form factor 
S(E) with the energy is unknown. However, as it is usually done in the literature (Fowler 
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et al. 1967, 1975; Wagoner et al. 1967), we can do a MacLaurin expansion: 

where ^| and are expressed in barn and barn Mev -1 respectively. The terms inside the 
brackets have no dimensions, therefore: 

5(0) ~ ^qcd (26) 

The dependence of the charged particle cross sections with a has been analyzed by Nollet 
and Lopez (2002), yielding: 

S(0) ~ a (27) 

Furthermore, it follows that all radiative capture rates should be multiplied by a factor t % ay , 
except the reactions T(cry) 7 Li and 3 He(o;7) 7 Be. This cross sections should be multiplied by 
f( a ) = Yl^i [ a toda V — l] (see table 2). Finally, in the cases in which the reaction produces 
two charged particles, the cross section should by multiplied by 1 — b + b at ? day (see table 3 ). 

We insert the expression for a(E) into the equation (21) in order to calculate the Boltz- 
mann averaged cross sections: 

where £ = 2rraZiZ 2 ^/ fic 2 /2kT and the masses in kg. The integrals are calculated in 
Bergstrom et al. (1999). Tables 4 and 5 show the dependence of charged particles reac- 
tion rates with the fundamental constants. 



3.1.1. Cutoff factor 

The truncated MacLaurin series we have use for S(E) diverges at high energy. Thus, 
it is important to include a cutoff factor for non-resonant reaction rates so that they can 



Table 2: Radiative captures, its dependence on a (f(a) = J2bi [ at od ay — l]) 



Reaction 


bo 


h 


b 2 


h 


h 


h 


3 H(a, 7 ) 7 Li 


1 


1.372 


0.502 


0.183 


0.269 


-0.218 


3 He(a, 7 ) 7 Be 


1 


2.148 


0.669 


-5.566 


-10.630 


-5.730 
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be used at any energy. The next term in the expansion for S(E) would be proportional to 
E 3 ~ T 2 , so as it is proposed in the literature (Fowler et al. 1967, 1975) we consider a cutoff 
factor: f co = e~^ T '^ Tco ^ 2 , where T co ~ and E r is the resonant energy (Fowler et al. 1967, 
1975). Therefore we multiply the expression (28) by a factor: 

fco = e"t * ) (29) 

This correction is relevant for the following reactions: 3 He(<i,p) 4 He, 3 H(rf, n) 4 He, 6 Li(p, a) 3 H, 
6 Li(a,p) 10 Be, 7 Li(p, a) 4 He. 



3.1.2. Alternative expression for the form factor 



In the MeV range the cross section form factor varies considerably. In this range the 
truncated MacLaurin series is not satisfactory so that it is convenient (Fowler et al. 1967, 
1975) to use for S(E) an expression of the form: S(E) = S(0)e~ aE . In such way, the cross 
sections are given by: 

(30) 



a (E) = c -aE ^olZ^J nc 2 12E 



E 

where a has no dependence on the fundamental constants. The quantities with units in 

equation 30 are 5(0) and E, therefore, in the chiral limit we have 5(0) ~ ^qcd- ^ n sucn 
way, the Boltzmann cross section (eq.(21)) yields: 



(av) 



8 / A y/2 S(0)kT (H 



^6 \hT) akT + 1 \A 



1/6 



,2N 1/3 

31 % 



1 + 



5 

36 



-1/3' 



cm" 



seg 



(31) 



Table 3: Dependence on a of different kinds of reactions rates 



Reaction 


Multiplied by 


Charged particles reaction rates 


a 

^today 


Photon emission 


a 

fytoday 


2 E(d,p) 


1 + 0.16-0.16^ 


3 He(n,p) 


1 - 0.30 + 0.30^ 


3 Re(d,p) 


1 + 0.09 - 0.09^ 


7 Li(p, a) 


1 + 0.18-0.18^ 


7 Be(n,p) 


1 - 0.20 + 0.20^ 
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where £ a = £y/akT + 1 = 2naZ 1 Z 2 y ^(akT + 1). This alternative expression for non- 
resonant reaction rates is relevant for the following reactions: 6 Li(p, 7) 7 Be, 3 H(a,7) 7 Li and 
3 He(a, 7 ) 7 Be 



3.2. Resonant charged particle reaction rates 

The expressions for the cross sections vary with the temperature. Moreover, in the range 
of energies relevant for our calculation there are certain reactions that proceed through many 
resonances. In this case, we have to include an extra term in the cross section. There are 
two kinds of resonances: i) Single Resonance, ii) Continuum Resonance. In each case we use 
the expressions given by Fowler et al. (1967, 1975). 



3.2.1. Resonance cross sections 

In this case, the following expression provides a good fit to the cross section: Fowler 
et al. (1975): 

, . 7rh 2 o; r r 1 r 2 

a(E) = 2^{E-E r )* + T*/A (32) 

where Tj is the partial with for the decay of the resonant state by the reemission of (i — 1) +i, 
T is the sum over all partial widths (the partial widths are not functions of a) , uj t = ( 1+5ab ' )9r 
and g r — 2 J r + 1, J r being the spin of the resonant state, \x in kg and E r is the resonance 
energy in the center of momentum system and depends on the nuclear radius. Finally, the 
Boltzmann cross section (av) can be calculated as follows: 

<™>=(^Y /2 ^* /fcT — (33) 
v 1 \ fikT J h seg K ' 

where 7 r = (^ 2 ) ■ Here, the cross section depends on the fundamental constants through 
the final temperature and the resonance energy. Besides, the resonance width is also a 
function of the fundamental constants, but the cross section is much less sensitive to this 
dependence. 

In a theory with massless quarks: E r ~ Aq C d- On the other hand, E r does not depend 
on a. The dependence of the temperature will be analyzed in section 4. This correction 
is relevant for the following reactions: 2 H(o;, 7) 6 Li, 6 Li(p, a) 3 H, 6 Li(o;, 7) 10 Be, 7 Be(p, 7) 8 Be 
and 7 Li(p, a) 4 He. 
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3.2.2. Continuum resonances 

When the temperature scale is of order Tg ~ 1 , there are several reactions that proceed 
through many resonances that are separated by intervals not greater than their widths or 
that overlap to form a continuum. In this cases, the cross section can be written as Fowler 
et al. (1975): 

t^o^-ir 1 ' 2 ^ e> c 

\ si E<C K J 

where m is integer or rational fraction, C is the effective continuum threshold energy and 
a(2C) is the cross section at E = 2C. 

After inserting this expression in the integral (21), we obtain: 



(av) = r(m + 3/2)a(2C)y^ ^V c / fcT (35) 

where T(m + 3/2) is the gamma function. On the other hand, C has units of energy. 
Therefore, in the chiral limit C ~ A.qcd- This correction is relevant for the following 
reactions: 3 He(<i,p) 4 He, 3 H(<i, n) 4 He, 6 Li(p, a) 3 H and 6 Li(a, 7) 10 Be. 



3.3. Non charged particles reaction rates 

In this case, there is no Coulomb barrier so the cross section cannot be written as the 
equation (23). Following Fowler et al. (1967) we write: 

a(E) = & (36) 
where v is the relative velocity. We consider the expression given by Fowler et al. (1967): 

s(B »^ (0)+ (I) fc / 1/2+1/2 (i)^ E (37) 

where E = E 1 ' 2 . 

In chiral limit o ~ ^qcdi an d therefore: 

5(0) ~ A Q 2 CD (38) 
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In this way, the expression for the reaction rates (equation 21) yields: 

<-^ s <°»( i+ ^S«^ i/2+ iS fcr ) (39) 

where ^§ is in units of s^MeV" 1/2 and in in units of s^MeV" 1 

dE se S dE se § 

For radiative emission reactions, the cross section should be multiplied by a factor —Sc^- 
Table 6 shows some reaction rates between a neutron and a nucleus. 

In some cases the reaction rates of the inverse reactions are needed. Next, we show the 
expressions for these reaction rates. For inverse reactions of the form [BCAn], where neither 
B or C are photons, we use the expression given by Fowler et al. (1967, 1975): 

[BCAn] = ^ + S ^)9A ( II^Y 2 e -Q/ k T [AnCB] J_ (4Q) 
(1 + dAn)9B9c \m B m c J seg 

where Q = tjia + rn n — mg — mc ■ In this case, the form factor 5(0) should be multiplied 
by a today because of the Coulomb barrier. 

For inverse reactions of the form [B7nA], we use the expression given by Fowler et al. 
(1967, 1975): 

UB.nA] = ( — r ( f§re-^W„)l- (41) 

(1 + S A n)9B V m B J V 2nh J se g 

where Q = m A + m n - m B , g n = 2j n + 1 = 2 and M v = 

In both cases, the additional dependence on the fundamental constants introduced by 
the inverse reactions proceed from the temperature (see section 4) and the masses (see section 
2). 



3.4. Neutron lifetime 

Neutron (3 decay is one of the few reactions whose cross section can be explicitly com- 
puted from first principles in terms of the fundamental constants. It can be approximated 
by the one point interaction of neutron, proton, electron and neutrino. The reaction rate for 
neutron (3 decay is: 

n -> v + e" + vl (42) 

Following Ichikawa and Kawasaki (2002) we write the inverse of neutron lifetime as 
follows: 

1 f Po 

d 3 p e d 3 p u S(E e + E u + m p -m n ) (43) 

r Jo 
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where, Gf is Fermi coupling constant, E e and E u are the electron and neutrino energy, and 
p e and p u are the electron and neutrino momenta. After integration we obtain: 

\ = D^ 7r2G ^ (v^ 3 ! (V - V - 8 ) + 15 9 ln (? + V 7 ?^)) ^ (44) 
where m e is the electron mass, q = — = m "~ T " p 



In such way, we obtain the dependence of the neutron decay rate r 1 on Gf and on the 
mass difference (which is a function of a and Aqcd, see section 2): 

4r = = 2—^ + 6.54 go - 3.839— (45) 



4. Abundances as functions of fundamental constants 

In this section we calculate the abundances of light elements and their dependence on 
fundamental constants. First we obtain the neutron abundance until the freeze-out time of 
weak interaction. After this time the neutrons decay freely into protons and electrons, so 
their abundance only changes due to this decay. 

The general form of the equations that govern the abundances of the light elements is: 

y. = j(t) _ V{t)Yl (46) 

where J(t) and T(t) are time-dependent source and sink terms and the dot corresponds to 
the time derivative. The time-dependent static solution of this equation is what we will call 
following Esmailzadeh et al. (1991) the quasi-static equilibrium (QSE) solution: 

/, = || (47) 

To determine the formation of light nuclei we shall solve the following equations using 
only the most important reactions according to the rates of production and destruction 
following the criteria established by Esmailzadeh et al. (1991): 

Y n = Y d Y d [ddn3] + Y d Y T [dTna] + Y p Y T [pTn3] + YdY^np] + 

-Y n Y p [npd 7 ] - Y n Y 3 [n3Tp] - Y n [n] (48) 

Y d = Y n Y p [npd 1 ]-2Y d Y d ([ddpT] + [ddn3})-Y d Y T [dTna} + 

-Y d Y 3 [d3pa] - Y d Yy [d^np] - Y d Y p [rfp3 7 ] (49) 

Y 3 = Y d Y p [pd3 1 ] + Y T Y p [pTn3} + Y d Y d [ddn3} + 
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-Y d Y 3 [dSpa] - Y n Y 3 [n3pT] (50) 
Y T = Y n Y 3 [n3 P T}+Y d Y d [dd P T}-Y d Y T [dTna} + 

-Y p Y T [pTn3] - Y p Y T [pT ia } (51) 

F 6 = Y d Y a [da^]-Y n Y6[n6aT\-Y p Y 6 \p6Ta] (52) 
Y 7 = Y n Y a [nQ7 1 ] + Y n Y B [nBp7}+Y T Y a [Ta7 1 } + 

-Y p Y 7 [p7aa] - Y n Y 7 [n78-f] (53) 
Y B = r p y 6 [p6S 7 ] + Y 3 Y a [3a£? 7 ] - Y^B-yZa] - y n y B [nSp7] + 

-Y P Y B [pB 7 8] - I^Fb [dBaap] (54) 

r Q = F d y 3 [d3po;] + Y n Y 3 [n3a^] + F d F T [dTna:] + Y p Y T [pT^a] (55) 

where n refers to neutron, p to proton, d to deuterium, T to tritium, 3 to 3 He, a to 4 He, 6 
to 6 Li, 7 to 7 Li, £? to 7 Be, 7 to the photon and [ijkl] is the rate of the reaction i + j — * k + l 
and Fj is the abundance of the % element relative to baryons (Yi = . In addition, these 
equations obey neutron number conservation: 

Y n + Y d + Y 3 + 2Y T + 2Y a = -Y n [n] (56) 



The method of Esmailzadeh et al. (1991) consists in calculating the different abundances 
between fixed point or stages. We shall solve equations 48 to 55 only for one element in 
each stage. For the other elements it is necessary to solve the quasi static equilibrium 
equation using only the most important rates of production and destruction. On the other 
hand, we perform the calculation of all final temperatures and abundances and all freeze-out 
temperatures numerically. Table 7 shows the different stages and the used equation. 

The equations that describe the production of n, D, 3 He and T are independent to the 
equations for 6 Li, 7 Li and 7 Be. Therefore, we shall solve the first three using the quasi static 
equilibrium equation and then we use these results to calculate the other abundances. 

To calculate the final abundance of light elements it is necessary to know the freeze-out 
temperature. The freeze-out of the production of each element happens when the most 
important destruction reaction rate equals to the expansion rate of the Universe. The 
dependence of the freeze-out temperatures and final temperature of each stage with the 
fundamental constants, will be calculated by deriving the equation that determines each 
temperature. 

Each section in this chapter will discuss the calculation of abundances during a certain 
stage. 
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4.1. Neutron abundance until the freeze-out of weak interaction, T > 9.1 x 10 9 K 



For the calculation of neutron abundance we follow the analysis performed by Bernstein 
et al. (1989). Let X pn (T) be the rate of weak process that convert protons into neutrons and 
X np (T) the rate of weak process that convert neutrons into protons. The basic rate equation 
reads: 



dX 

— = \ pn (t){l-X(t))-\ np (t)X(t) 



(57) 



where t is the time, and X is the ratio of the number of neutrons to the total number of 
baryons. After changing variables (y = ^p) , the solution of the last equation can be written 
as follows: 



X(y) = X eq (y) + dyey X eq (y ) 
Jo 



D K(y)-K{y ) 



(58) 



where 



K(y) = b 
X eq (y) ~- 



^ + ^7 + -+ I ^ + ^7 }e~ y 



v y 



y 3 y 2 



b = 255 



P i 



45 



Am 2 r V 43tt 3 



1 + ev 



(59) 



r is the neutron mean life and Am = m n — m p . In order to obtain the asymptotic behavior, 
the limit T — > or y — > oo is taken: 



X(y = oo) 



dy e y X eq [y ) e 



-K 



(0 



0.151 



(60) 



In the last equation, only b depends on the fundamental constants through r and Am 
(see sections 2 and 3 for the dependence of these quantities with the fundamental constants). 
In such way, from equation 60, we obtain: 



5X{y = oo) 
X(y = oo) 



= -1.04 



6G, 



-2.361^^ + 1.386^ 



5a 



G, 



QCD 



(61) 



4.2. Until the production of 4 He becomes efficient, 9.1 x 10 9 K > T > 0.93 x 10 9 K 

After the freeze-out of the weak interactions, the only change in the neutron abundance 
is due to neutron decay. Therefore, the neutron abundance in this stage reads: 



Y n = X(y = oo) e-^ = e' ' 198 ^ 



(62) 
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In the beginning of this stage there are no nucleus with two or more nucleons, therefore 
it is a good approximation to consider: Y p ~ 1 — Y n . However, as the universe expands, 
the temperature goes down and light nuclei formation begins. Therefore, at the end of this 
stage, this expression is no longer valid. 

In order to get a consistent solution of equation 56 (Esmailzadeh et al. 1991), it it is 
necessary to set all the rates equal to zero with the exception of the largest rate which equals 
to — 2Y a — Y n [n\. In such way, the equations to solve in this stage are: 

Y n = -2Y a -Y n [n] (63) 
Y d = Y 3 = Y T = (64) 

Table 8 shows the solutions. 

When the production of 4 He becomes efficient the stage ends. The final temperature is 
given by is given by setting Y n = in equation 63. For this stage, we obtain T/ = 0.93 and 
the following results: 

Y p f = 0.76 Y d f = 4.1 x 10" 4 Yf = 2.0 x 10~ 5 

Y r { = 0.12 Y 3 f = 5.8 x 10~ 8 y£ = 0.06 

where Y/ is the final abundance of each nucleus or nucleons on this stage. It follows that 
the the abundances of D, T and 3 He are negligible respect to the abundances of neutrons 
and 4 He. This means: 

Yj = \-Yl-Yl-Yi-1Yl-7Yl 

~ \-Yl-TYl = \-TY,l (65) 

Now, in order to calculate the dependence of the final temperature with the fundamental 
constants for this stage, we derivate the equation 2Y a = Y n [n] with respect to the fundamental 
constants and the temperature. In such way, we obtain: 

5Tl 5tt B h 2 5G F 5a 5A cd , N 

— i. = 0.068 -0.053— ^ + 0.063— + 0.871 - QC 66 
T( Sl B h 2 G F a Aqcd 



where we also considered the dependence with the baryon fraction. Finally, the dependence 
of the abundance of neutrons on the fundamental constants and fish 2 yields: 

^ = 0.029^ - 1-522^ + 2.296^ - 3.459^ (67) 
Y n tt B h 2 G F a Aqcd 
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4.3. Final abundance of 4 He 

The next stage corresponds to the calculation until the rate of production of deuterium 
dominates over the rate of change of neutrons. However, the freeze-out temperature of 4 He 
(T = 0.915 x 10 9 K) is lower than the final temperature of the previous stage but bigger than 
the final temperature of the next one. Therefore, we calculate now the final abundance of 
4 He. In this case, the neutron number conservation equation reads: 

2Y a = Y n [n\ (68) 

For the others nucleus the quasi static equilibrium equation is valid (see table 8). The 
production of 4 He is dominated by [dTna] and [pTja] : 

Y a = Y d Y T [dTna] + Y^pT^} = (^ /[^j?] ) [MpT] (69) 

After solving numerically for Tg the equation (68), we obtain Tg = 0.915 and Y[ = 
2Y n = 0.238. When the rate of 4 He production equals to the rate of neutron destruction, 
there is no more neutron that can form 4 He. Since this happens earlier than the usual 
freeze-out-time, we use equation 68 to calculate the freeze-out temperature. In such way, 
the dependence of the freeze-out temperature on the fundamental constants and Qsh 2 yields: 

^ = 0-061^C " 0-052^ + 0.063^ + 0.869^ (70) 
T 9 a tt B h 2 G F a A Q cd 1 ' 

Finally, since Y£ = 2Y n , we can express the variation of the final abundance of 4 He as 
a function of fundamental constants and fi^/i 2 : 

^ = 0.029^ - 1.538^ + 2.324^ - 3.496^ (71) 
Y£ Q B h 2 G F a Aqcd 



4.4. Neutron cooking, 0.93 x 10 9 K > T > 0.765 x 10 9 K 

In this section we shall calculate the deuterium abundance as long as the change of 
neutron dominates the deuterium production rate. This is valid until the production rate of 
deuterium dominates the rate of change of neutrons, so this stage is over when Y n = Y d . In 
this stage, the neutron number conservation equation reads: 



Y — —2Y 

1 n — ^ 1 a 



(72) 
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For D, T and 3 He we solve the quasi-static equilibrium equations. The solutions are shown in 
table 8. For 4 He we solve the complete equation but considering only the largest production 
term Y d Y T [dTna\. Inserting all these solutions in equation 72, we obtain: 

Y n = M Y n Y p ^L) 2 [ddpT] (73) 



"Y^np]) 

where the initial condition is given by the final values of the previous stage: F„° = 0.12 and 
T 9 ° = 0.93. We can write the solution to the last equation as follows: 

After changing the integration variable to T 9 we perform the integral numerically as a func- 
tion of temperature. We also compute the final temperature of this stage using the condition: 

Y n = Y d (75) 

We obtain: 

T/ = 0.765 Y n = 6.4 x 10" 4 = Y d (76) 

From 75 we obtain the dependence of the final temperature of this stage with respect to the 
fundamental constants and fish 2 : 

H = 0-031^# + 0.015^ - 0.023^ + 1.034^ (77) 



T/ ' V B h 2 G F a A 



QCD 



Finally, the dependence of the final neutron and deuterium abundance can be obtained from 
equation (74): 

5Yd 6Yn -1.099^- 0.058^ + 1.871^-0.488^ (78) 



Y d Y n ' n B h 2 G F a A. 



QCD 



4.5. Deuterium cooking, T — > 

For temperatures lower than T 9 = 0.765, the largest production rate corresponds to 
deuterium. Therefore, we set all other derivatives to zero in equation 56. Since the largest 
term for deuterium destruction is tritium production, the equation to solve is: 

Y d = -2Y d Y d [ddpT] (79) 
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with the initial condition F d ° = 6.4 x 10 4 on T 9 ° = 0.765. Since this equation has the same 
form of equation 73, the solution reads: 

Y d=(yo+ 2 [ [dd P T]dt) (80) 

V d J t initial / 

In order to calculate the deuterium final abundance we consider the limit T — > ( 
t — > oo). We obtain numerically, the deuterium final abundance Yj = 2.410 x 10~ 5 . On 
the other hand, the dependence of the deuterium final abundance with the fundamental 
constants and fish 2 can be calculated by deriving equation 80: 

5Y? 5tt B h 2 5G F 5a SA OCD , , 

= -1.072— B — - 0.036— ^- + 2.320— + 0.596 QCD 81 
Yd U B h 2 G F a Aqod 



4.6. Final abundances 

Here we calculate the freeze-out temperature and final abundances of 3 He, T, 6 Li, 7 Be 
and 7 Li and the dependence of these quantities with the fundamental constants. In order 
to calculate any light element abundance it is necessary to solve the quasi-static equilibrium 
equation: 

Y t = (82) 

We solve these equations considering only the most relevant reactions. In table 9 we show 
the quasi-static equilibrium solutions. 

In order to compute the freeze-out temperature, we set the largest rate of destruction 
T of each equation that governs the abundance of the light elements equal to the universe 
expansion rate H: 

T = H= T^Tlseg- 1 (83) 

Table 10 shows the different freeze-out temperatures and their dependence on fundamental 
constants which is calculated deriving the previous equation. Using the freeze-out tempera- 
ture we calculate the final abundance of the different nucleus and their dependence on the 
fundamental constants and VLsh 2 . In table 11 we show these results. 



5. Results and discussion 



In this section we compare the theoretical predictions of the abundances of the light 
elements obtained in the last section with observational data. 
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In section 4 we have obtained 7 equations of the form: 




+ A 



5Q B h 2 
Q B h 2 



(84) 



where Aj, Ci and Dj are constant coefficients (see table 11) 



' Y 
1 1 



— ySBBN 



and 



Y.sbbn anc | Y° bs are the theoretical and observed abundance respectively. 

However, independent observational data are only available for the abundances of D, 
3 He, 4 He and 7 Li. In table 13 we show the independent data we consider in this work. For a 
recent review of all observational available data on primordial abundances see Particle Data 
Group et al. (2004). On the other hand, recent papers (Coc et al. 2004a; Cyburt 2004) have 
brought the attention to the errors introduced by the values of the cross sections involved in 
the calculation of the abundances. Cyburt (2004) has also analyzed the propagation through 
the theoretical abundances, yielding a "theoretical" percent error of 5%. In the original work 
of Esmailzadeh et al. (1991), the error introduced by the semi-analytical method is estimated 
to be of order 5%. Therefore, we will add in order to solve system 84, to the errors of table 
13 an error of order 10%. 

First we perform a test to check the consistency of the data (Riveros and Vucetich 
1986). For each group of data (Y,j) belonging to the same abundance, we calculate the 
weighted averaged value Y and its corresponding error <jj. Then we compute: 



If the errors are Gaussian distributed, the expected value of x 2 is (k — 1) where k is the 
number of data in each group. Furthermore, the corresponding ideogram of each group of 
data (Particle Data Group 2002), should be a Gaussian. It follows from figure 1 and from 
the calculation of x 2 that D and 4 He data are not Gaussian distributed. However, since 



6 = y ^-j- is not that greater than one, we can use the data but increasing the observational 
error by a factor 0. The values of 6 are 2.4 for D, 2.33 for 4 He. 

We assume that any difference between the theoretical abundance and the observational 
abundance is due to the variation of fundamental constants. In such way, the solution of 
system (84) gives a constraint to this variation. The solution is given by (Arley and Buch 
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1968): 




(86) 
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where B is the n x 4 matrix, n is the number of observational data: 



B = 



S is the n x 1 matrix: 



/ ^ fix d £>! \ 

Ai Bo Co Do 



\ -Bn Cn -D« / 



V 2 



(87) 



\ §Xn I 



and P is the n x n matrix of weight: 



f Pi ... \ 
p 2 ... 



\ 



(89) 



Pn J 



where pi 



\ and <7i are the observational errors. 



The most accurate estimation of Qsh 2 arrives from constraining parameters with data 
from the CMB provided by WMAP (Spergel et al. 2003). Fixing the baryon fraction with the 
WMAP value (i.e. setting Dj = 0), the results of solving the system 84 with all data listed 
in table 13 are shown in table 14. These results are consistent within 3a with variation of 
the fundamental constants. On the other hand, the results considering only variation of the 
fine structure constant are shown in table 15. These results are consistent with no variation 
of a within 3c. In order to rule out any systematic error of the data, we computed the 
solution of system 84 again but excluding one group of data at each time. Again, the results 
are consistent with variation of the fundamental constants in all cases but in the case where 
the 7 Li data were excluded (see tables 14 and 15). 

Even though, the WMAP estimate of the baryon density is the most accurate one, it 
is still affected by degeneracies with other cosmological parameters (Spergel et al. 2003). 
Therefore, we added an independent estimation of f2#/i 2 in our analysis. In appendix A we 
use data from X-ray measurements, galaxy surveys and cepheids calibration in order to get 
an independent value of the baryon density. Furthermore, we computed again the results 
of sections 3 and 4, changing the value of Vl B h 2 to 0.0223. This value is the weighed mean 
value between the WMAP estimate and the value of appendix A. However, we found no 
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difference in the value of the coefficients of the variation of fundamental constants and Vt B h 2 . 
The results obtained solving system 84 including both estimates for the baryon fraction (i. 
e. Di ^ ) show no significant difference with respect to the case where only the WMAP 
value was considered (see tables 16 and 17). Furthermore, in order to check for consistency 
of our method, we solved again system 84 allowing only for variation of fish 2 with respect 
to the weighed mean value (i.e. Ai = Bi = Gi = 0). These results are shown in table 18. 
On the other hand, in order to learn about the degeneracies of the fundamental constants 
within the BBN model, we computed the correlation coefficients from the error matrix. We 
find that there is high correlation between a and Aqcd, ol and Gf and Aqcd and Gf, while 
the correlation between other pairs of parameters is not significant. 

In order to understand the discrepancy of the results obtained with and without the 
7 Li data, we computed the relative residuals (Arley and Buch 1968), and their respective 
theoretical and empirical probability in both cases. Figure 2 shows that in the case where 
both the variation of the fundamental constants and the deviation of Vtsti 2 from the WMAP 
estimate is considered, the theoretical and empirical probability distributions are very similar, 
while in the case where only the deviation of VLbIi 2 is considered, there is slight difference 
between the empirical probabilities (both with all data and excluding 7 Li data) and the 
theoretical probability. Including the variation of fundamental constants gives more degrees 
of freedom to system 84. Therefore, we suspect that the possible non reported systematic 
uncertainties "hide" under the variation of the fundamental constants. On the other hand, 
we performed a Kolmogorov-Smirnov (K-S) test, in order to check the goodness of our fit. For 
the results obtained considering variation of all constants and fish 2 , we obtain a probability 
of 21% to obtain a worse fit, while excluding the 7 Li data the probability lowers to 11%. 
On the other hand, if we only consider the deviation of Vtsh 2 with respect to the WMAP 
data, we obtain a probability of 99% for all data, while excluding the 7 Li data gives a 49% 
of probability to get a worse fit. However, we consider the results of the K-S test only 
indicative, since even though the data considered are independent the residuals are not. 

We mentioned in the introduction the disagreement between the 7 Li observational abun- 
dances with the D observational abundance and WMAP estimate of the baryon density. 
Richard et al. (2005) claim that a better understanding of turbulent transport in the stars 
is necessary to understand this discrepancy. Moreover, Melendez and Ramirez (2004) have 
reanalyzed the 7 Li data with an improved infrared flux method temperature scale, obtaining 
values that are marginally consistent with the WMAP estimate. However, solving system 
84 with the 7 Li abundance taken from their work, does not change in a significant way our 
results. 

We adopt the conservative criterion that the third and fourth column of tables 14 and 



-25 - 



16 are the constraints on the variation of the constants we obtain with the method and 
hypothesis described in this paper. We also consider that more observations of 7 Li are needed 
in order to arrive to stronger conclusions. However, if the all data are correct, this analysis 
shows that varying coupling constants may solve the concordance problem between BBN 
and CMB. Our results within 2 a are consistent with the analysis performed by Ichikawa 
and Kawasaki (2004), where only the variation of a and a non standard expansion rate. 



A. Appendix I 

In this appendix, we combine independent astronomical data in order to obtain and 
independent estimation of the baryon density. From measurements of hot gas in clusters it 
possible to obtain an estimate of ^-/i 3 / 2 . 

Ettori (2003) has brought the attention to the fact that the contribution from baryons 
in galaxies and "exotic sources" like intergalactic stars and baryonic dark matter are not 
considered in the results obtained from measurements of hot gas in clusters. Furthermore, 
Donahue et al. (2003) have estimated the contribution from the galaxies as follows: f ga i = 
0.15/i 3//2 fgas while the "exotic" contribution has been estimated in f eX otic = 0.3 f ga i (Ettori 
2003). Therefore, we add to the estimation of the baryon fraction done by Donahue et al. 
(2003) the contribution from galaxies, yielding the following value: 

9lLh 3/2 = 0.0737 ±0.0143 (Al) 

^ 'm 

The values of the other estimates (Majerowicz et al. 2002; Castillo-Morales and Schindler 
2003) are contained within the error in this estimation. 

On the other hand, Q m h has been estimated from large redshift galaxy surveys like 
Sloan Digital Sky Survey (Pope et al. 2004) and 2dF Galaxy Redshift Survey (Percival 
et al. 2001), while the most stringent bound on the Hubble constant follows from cepheid 
calibration (Freedman et al. 2001). Thus, combining all these data (see table 19) , and after 
propagating errors, we obtain the following value for the baryon density: 

n B h 2 = 0.017 ±0.007 (A2) 

This value is less accurate that the estimation done with the data of WMAP (Spergel 
et al. 2003) 

tt B h 2 = 0.0224 ± 0.0009 (A3) 
but we will consider it in order to have an independent data of this quantity. 
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Table 4: Charged particles reaction rates G = f2£/i 2 T 9 a 1 / 3 /! 1//3 , ^ = /ia 2 , S(6) 
tt B h 2 fi- b , S(a) = (l + a-a^^), P /B (#, T 9 , ci, c 2 , c 3 , c 4 , c 5 ) = 1 + Cl 

1Q -12 ^-1/3 T V3 + ^ x 10 10^1/3 T 2/3 + CsTg + Q x l 20 ^2/3 T 4/3 + ^ x 1Q 10 ^1/3^/3 



Reaction 


Reaction rate 






3 H(p, 7 ) 4 He 


1.14 xHT 7 [S(0)] 2 e UJ x 
P IB (^, T 9 , 4.36, 4.14, 1.26, 3.35, 2.61) 


2 H(d,n) 3 He 


2.26 x 10~ 3 E(0) e " 9 - 55xl ° 10 (^) x 
P IB (V, T 9 , 4.36, 1.96, 0.6, -0.206, -0.16) 


2 H( 3 He,p) 4 He 


0.39 £(0.09) e - L52x10 (^) -^ T ^\ 

Pib(V, T 9 , 2.75, -2.16, -0.42, 13.5, 6.58) + 
+2.63 x 10~ 8 S (|) £(0.09) T 9 5/2 e" 1 - 76 ^ 1 


3 H(d,n) 4 He 


0.49 £(0) e -^xio-(^) 1/3 -(ii4i.67T 9Q ) 2>< 

P IB (^, T 9 , 4.36, 3.78, 1.16, 46.8, 3.64 x 10 11 ) + 
+3.39 x 10~ 8 S (|) £(0) T 9 7/3 e" - 523 ^ 1 


2 H(d,p) 3 H 


2.37 x KT 3 £(0.16) e - 9 - 55xl0l °(^) i/3 x 
P IB ($>, T 9 , 4.36, 1.16, 0.35, -0.051, -0.04) 


2 H(«, 7 ) 6 Li 


1.88 x 10- 10 [£(0)f e- 1 - 5 ^ 101 ^^) 1 '^ 

P IB (V, T 9 , 2.75, -9.9, 8.85, -2.43, -1.19) + 
+8.27 x 10~ 39 S (|) [£(0)f T 3/2 e" 8 - 2287 ^ 


H( 6 Li,a) 3 H 


0.20 £(0) e -^xio"(f) 1/a -(24.94T 9Q ) 2>< 

P/b(^, ?9, 2.10, -0.14, -0.02, 0.033, 0.012) + 
+4.53 x 10~ 8 S (|) £(0) T 2 e - 21 - 82T 9 _1 + 
+6.68 x 10~ 34 S (|) E(0) T 3/2 e" 17 - 76 ^ 1 
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Table 5: Charged particles reaction rates = Qb^Tq 3 a 1 / 3 fj, 1//3 , \1/ = /ia 2 , = 
tt B h 2 fi- b , S(a) = (l + a - a-^) , P /B (^, T 9 , ci, c 2 , c 3 , c 4 , c 5 ) = 1 + c x x 

10 -12^-l/3 T V3 + ^ x 10 10^1/3 T 2/3 + C3Tg + Q x l 20^2/3 T 4/3 + ^ x 1() 10 ^1/3 ^ 

P LN (x,d 1 ,d 2 ,d 3 ,d 4 ,d 5 ) = 1 + rfii + rf 2 a; 2 + d 3 :r 3 + rf 4 x 4 + rf 5 x 5 , T 9t , = vT 9 + 11.605 



Reaction 


Reaction rate ( — ) 


3 H(a, 7 ) 7 Li 


[E(0)] 2 P LN (^?, 3.17, 0.50, 0.18, 0.27, -0.22) x 
|7.47xKT 6 e - L52xl0ll (^) 1/3 x 

P IB (*, T 9 , 2.75, -0.76, -0.15, 0.36, 0.18) + 

+2.68 XlO" 5 T -5/6 e -6.69xl0-(^) | 

v = 1.59 


3 He(a, 7 ) 7 Be 




{ 

H 

f 


P(0)] 2 P LN (^^, 2.15, 0.67, -5.57, 10.63, -5.73) x 
's.27 xlO- 5 e" 2 - 41xl0ll (^) 1/3 x 

P IB O, T 9 , 1.73, -0.0019, -0.00024, -0.00028, -8.8 x 10~ 5 ) + 
-3.12xl0-X 5/8 e" 1 - 06xl0ll ("2") / } 
= 1.24 


H( 7 Li,a) 4 He 


3.33 xKT 3 S(0) e -^x^(f) 1/a -(o. 22 ^ x 
Pib T 9 , 2.10, 3.65, 0.54, -5.30, -1.98) + 

+2(1) E(0) T 9 3 / 2 x 
[5.54 x io- 34 e- 30 - 44 / T 9 + 7.98 x lo-^e" 4 - 479 /^] 


H(d, 7 ) 3 He 


1.11 xlO" 8 [S(0)] 2 e - 9 - 545xl0l0 (^) 1/3 x 
Pib T 9 , 4.36, 8.66, 2.65, 1.26, 0.98) 
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Table 6: Non charged particles reactions rates, S nc (a) = 1 + a — a at " day , Pf(T9,Ci,c 2 ) = 

-1 . Til/2 1 m m n —m v 

l + CiV +c 2 T 9 , g= p 



Reaction 


Reaction rate ( — ) 


n( ,e")H 


§n 2 G 2 F m 5 e [Vg 2 - 1 (2g 4 - 9q 2 - 8) + 
+ 15g In (q + y/q 2 - l) 


H(n, 7 )d 


40.92 ft B /i 2 S„ c (-1) T 9 3 P F (T 9 , -0.86, 0.43) 


F 7 {H( 7 ,n)H} 


2.70 x 10 49 S nc ( 1) KM 3/2 T 9 3/2 e -n-605 ed /T 9>< 
P F (T 9 , -0.86, 0.43) 


3 He(n,p) 3 H 


6.53 x 10 5 tt B h 2 T 3 S nc (0.3) P F (T 9 , -0.59, 0.1832) 


3 H (p,n) 3 He 


/ \ 3 / 2 

6.53 x 10 5 Q B h 2 S nc (0.3) E nc (-1) (»J x 
Pf(T 9 , -0.59, 0.1832) e - n - 605 ^/^T 9 3 


7 Be (n,p) Y Li 


6.27 x 10 6 Vl B h 2 S nc (0.2) T 3 P F (T 9 , -0.903, 0.215) 



Table 7: Stages and equations 



Stage 


Equations 


Final temperature 


Until the weak interaction freeze-out 






Until the production of 4 Hebecomes efficient 


Y n = -2Y a -Y n [n] 
Y d = Y 3 = Y T = 


2Y a ~ Y n [n] 


Production of deuterium dominates the 
rate of change of neutrons 


Y — —TV 
Y d = Y 3 = Y T = 


Y n = Y d 


Deuterium final abundance 


Yd = — 2F Q 
Y n = Y 3 = Y T = 


T 9 ^0 



Table 8: Solutions of the quasi static equilibrium equations for each stage 



T 9 


Nucleus 


Solution 


9.1 > T 9 > 0.93 


D 


Y — v V ^™ pd7 ' 

*d 1 n 1 p Y [d-ynp] 


T 


Y y d y P \pd-dy\ + v d y d [ddn3\ + y T v p \ P >i W3\ 


lT Y d \d3pa\+Y n \n3pT] 


3 He 


Y - y n y s \n3pT\+y d y d \ddpT\ 

3 Y d [dTna] +Y P [pTn3] 


0.93 > T 9 > 0.765 


D 


Y — V V \- n v d "i\ 

id 1 n 1 VY~,\d~ i np\ 


T 


Y 3 = Y d - 


ddn3 
d3pa 




3 He 


Yt = Y d - 


ddp'l 
dTna 


J 
1 


T 9 ^0 


D 


Y — V V i n P d l\ 

ld InI PY y [d-ynp] 


T 


Y 3 = Y d - 


ddrd 
d3pa 




3 He 


Y T = Y d - 


ddp'l 
dTna 


J 
1 



Table 9: Quasi-static equilibrium solutions 



Nucleus 


Quasi-static equilibrium solutions 


3 He 


Y _ Y d \ddn3\ 
3 [d3pa\ 


T 


Y y d WpT\ 


6 Li 


Y _ Y d Y a [da6j\ 
6 Y p \p6Ta] 


7 Be 


Y Y 3 Y a [3aBj\ 
IB Y n [nBp7] 


7 Li 


Y Y n Y B [nBp7\ + Y T Y a [Ta7j\ 
17 Y p [p7aa] 



Table 10: Freeze-Out temperature and their dependence on fundamental constants, 

w 5Gf + R Sa + T S AQ£b + jSUsh? 



Nucleus 


Equation 


rpFreeze—Out 
1 9 


W 


R 


T 


J 


3 He 


Y d [d3pa] = H 


0.403 


0.008 


-0.510 


1.168 


0.016 


T 


Y d [dTna] = H 


0.105 


0.009 


0.122 


1.181 


0.018 


6 Li 


Y p [ P 63a] = H 


0.069 


-0.076 


1.962 


1.118 


-0.156 


7 Be 


Y n [nBp7] = H 


0.319 


0.217 


-0.712 


1.39 


0.350 


7 Li 


Y p \p7aa] = H 


0.185 


-0.088 


1.692 


0.946 


-0.182 
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Table 11: Abundances and their dependence on fundamental constants, = A-^ + B^- + 

Aqcd n B h 2 





Abundance 


A 


B 


C 


D 


2 H 


2.741 x 1(T 5 


-0.036 


2.320 


0.596 


-1.072 


3 He 


6.95 x 1(T 6 


-0.051 


0.983 


0.999 


-1.102 


3 H 


1.21 x 10" 7 


-0.041 


0.252 


0.941 


-1.083 


4 He 


0.238 


-1.538 


2.323 


-3.497 


0.029 


6 Li 


5.7 x 10" 14 


-2.061 


7.414 


-3.462 


-1.047 


7 Be 


5.60 x 10~ 10 


-0.172 


-9.450 


-1.038 


2.209 


7 Li 


2.36 x 10~ 10 


-0.720 


1.824 


-3.411 


0.068 



Table 12: Theoretical abundances in the standard model the WMAP estimate Vl B h 2 = 0.0224 



Nucleus 


ySBBN _|_ fiy. SBBN 


2 H 


(2.51 ±0.37) X 10~ 5 


3 He 


(1.05 ±0.15) x 10~ 5 


4 He 


0.2483 ±0.0012 


7 Li 


(5.0 ±0.3) x 10- 10 
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Table 13: Observational abundances used in this work 



Nucleus 


Yobs _|_ fiY obs 


Cite 


D 


(1.65 ±0.35) x 10~ 5 


Pettini and Bowen (2001) 


D 


(2.54 ±0.23) x IO" 5 


O'Meara et al. (2001) 


D 


(2.42lg;») x IO" 5 


Kirkman et al. (2003) 


D 


(3.25 ±0.3) x 10~ 5 


Buries and Tytler (1998a) 


D 


(3.98+ , 5 !) x IO" 5 


Buries and Tytler (1998b) 


D 


(l.6lg;») x 10-5 


Crighton et al. (2004) 


3 He 


(1.1 ±0.2) x IO" 5 


Bania et al. (2002) 


4 He 


0.244 ±0.002 


Izotov and Thuan (1998) 


4 He 


0.243 ± 0.003 


Izotov et al. (1997) 


4 He 


0.2345 ±0.0026 


Peimbert et al. (2000) 


4 He 


0.232 ±0.003 


Olive and Steigman (1995) 


7 Li 


(1.23+0-1) x IO" 10 


Ryan et al. (2000) 


7 Li 


(i-S8±8:lS) x io- 10 


Bonifacio et al. (1997) 


7 Li 


(1.73 ±0.05) x IO" 10 


Bonifacio and Molaro (1997) 


7 Li 


(2.191^) x IO" 10 


Bonifacio et al (2002) 




-4-3 -2 -1 1 2 3 4 5 -4-3 -2-10 1 2 3 4 



Fig. 2. — The full line shows the theoretical probability of the residuals, the dotted line shows 
the empirical probability computed with all data and the dotted line shows the empirical 
probability computed with all data but 7 Li. Left: Only deviation of VLsh 2 with respect to 
its mean value is considered; Right: variation of all constants and deviation of f^/i 2 from 
the mean value is considered. 



-40- 



Table 14: Constraints on the variation of fundamental constants (Vl B h 2 = 0.0224). 





All data 


All data but 7 Li 




Value 


a 


Value 


a 


SG F 
G F 


-0.886 


0.053 


-0.257 


0.659 


Sa 
a 


-0.136 


0.041 


-0.054 


0.097 


dA QCD 
A QCD 


0.309 


0.023 


0.087 


0.233 



Table 15: Constraints on on the variation of a (Vl B h 2 = 0.0224). 





All data 


All data but 7 Li 




Valor 


a 


Value 


a 


Sa 
a 


-0.041 


0.024 


-0.015 


0.005 



Table 16: Constraints on the variation of the fundamental constants using two independent 
estimates for the baryon fraction. 





All data 


All data but 7 Li 


Value 


a 


Value 


a 


8£l B h z 

n R h 2 


0.004 


0.036 


0.0005 


0.039 


6G F 

G F 


-0.886 


0.050 


-0.258 


0.64 


Set 
a 


-0.134 


0.044 


-0.053 


0.095 


S^QCD 


0.310 


0.023 


0.087 


0.229 



Table 17: Constraints on the variation of a using two independent estimates for the baryon 
fraction. 





All data 


All data but 7 Li 


Value 


a 


Value 


a 


Sa 
a 


-0.086 


0.034 


-0.015 


0.005 
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Table 18: Constraints on the deviations of Qsh 2 respect to the mean value considered in this 
work (0.0223). 





All data 


All data but 7 Li 


Value 


a 


Value 


a 


sn B h 2 
n B h 2 


-0.085 


0.294 


-0.014 


0.054 



Table 19: Observational data used to perform an estimate of the baryon density 



ns h 3/2 


Cite 


0.067 ±0.03 
0.073 ±0.013 
0.056 ± 0.007 


(Donahue et al. 2003) 
(Majerowicz et al. 2002) 
(Castillo-Morales and Schindler 2003) 




Cite 


0.20 ±0.03 
0.207 ±0.030 


(Percival et al. 2001) 
(Pope et al. 2004) 


h 


Cite 


0.72 ±0.08 


(Freedman et al. 2001) 



